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Abstract 

Based on the t — J model and the self-consistent Born approximation, the 
damping of quasiparticle hole states near the Fermi surface is calculated in 
a low doping regime. Renormalization of spin-wave excitations due to hole 
doping is taken into account. The damping is shown to be described by a fa- 
miliar form ImS(k',e) oc (e'^/ei?) ln(e/ei?) characteristic of the 2-dimensional 
Fermi liquid, in contrast with the earlier statement reported by Li and Gong 
[Phys. Rev. B 51, 6343 (1995)] on the marginal Fermi liquid behavior of 
quasiholes. 
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Understanding of the quasiparticle (QP) characteristics of charge carriers forming the 
normal-state electronic properties of high-T^ superconductors (HTSC's) is an issue of current 
interest .0 In particular, it is now under debate whether these compounds can be described 
within the normal Fermi liquid (FL) approach, or a more exotic scenario, for instance, the 
marg inal FL (MFL) concept,! should be involved. 

In an attempt to understand the QP properties of HTSC's one has to take into ac- 
count a strong difference between an intermediate and a low doping regime as it follows 
from the angle resolved photoemission (ARPES) experiments.i Actually, at the intermedi- 
ate level of doping ARPES indicates a large Fermi surface (FS), while the reference insulator 
compoundd^ show a hole dispersion that is compatible with a small four-pocket shape of FS 
at low doping. In the latter regime, a hole propagation is strongly affected by the presence 
of antiferromagnetic (AFM) correlations in the spin subsystem. The essential features of 
this problem are described by the t — J model. 

In the present paper, we investigate the regime of low doping based on the t — J model 
in the slave-fermion Schwinger boson representation. We are mainly interested in a QP 
hole behavior near FS. We show that at zero temperature the imaginary part of the hole 
self-energy ImE(k, cj) oc cu^lncj, which indicates a conventional FL behavior of quasiholes. 
Our result is at variance with one reported in Ref. ^, where the MFL-behavior of quasiholes 
is obtained even at T = 0. The reason of this contradiction is discussed below. 

By using the slave-fermion Schwinger boson factorization for electron operators, the t — J 
model with the Neel ground state can be mapped onto the so-called spin-polaron§^ model 
with the Hamiltonian given by 



q 



+4= E 4^k-q[Mi(k, q)«q + M2(k, q)«lq], (1) 

V A k,q 



where cUq = — 7q with 7k = \{coskx + cosfc^), /i is the chemical potential of holes, 

Mi(k, q) = M2(k — q, — q) = ('Uq7k-q + fq7k), and z = 4 for the square lattice; A^ is the 
total number of sites, the lattice spacing is taken to be unity, and Wq, fq are the usual 



parameters of the Bogoliubov u — v transformation. In the Hamiltonian (|l|) /ik {h\^ are 
canonical spinless fermion operators and ctq (a^) are canonical boson operators. 

We introduce the Fourier transformed two-time retarded Green function (GF) 
G(k, uo) =<^ /ikl^k for fermions and the matrix GF -D(q, u) =^ ^ql^q for magnons, 
y4q is the two-component operator and A]^ = (a^, a-q). By applying the irreducible GF 
method^ and using a decoupling procedure, which is equivalent to the self-consistent Born 
approximation (SCBA), both for G(k,uj) and D(q, cj), we obtain 

G-\k,Lo) = uj + n-E{k,uj), (2) 

uj^ + u + U22{q,uj) -ni2(q,u;) 

-1121 (q,t^) t^q - u; + nii(q,a;)^ 

where Pq(u;) = [cu - n"(q,cj)]2 - [cjq + n+(q,cj)]2 + ni2(q, u;)n2i(q, tu) with n=^(q,u;) = 
l/2[nii(q, cij) ±n22(q, uj)]- The elements of the polarization operator n(q, uj) for the magnon 
GF has the form 

(zt) ff 

na/3(q,w) = ^^X!^"/3(k,q) jj duidu^ 

^ — oo 

X - n{uj{)] ' [ , 3 

- tUi + CJ2 + 

where 5fa^(k, q) = MQ,(k, q)M^(k, q) and p\<jyuj) = — l/7rImG(k, cu) is the hole spectral func- 
tion. The hole self-energy S(k, cu) is given by 

^ — oo 

^^ Pk-q(^l)Xk,q(^2) /^v 

to" — cji — 0^2 + ir] ' 

where n{ijj) = (e^'^ + 1)^^ and A^(ti;) = (e^'^ — l)^"*^. In (Q) we have introduced an effective 
spectral function Xk,q(^^) as 

Xk,q(^) = --X!^a/3(k,q)ImL'a/3(q,^), (5) 

a/3 



for spin fluctuations coupled to a particular hole k-state. 

For the single-hole problem at zero temperaturellll the self-energy (^) is reduced and only 
the contribution due to Dii{ci,uj) remains. Considering the case of finite hole concentration 
6, the full form (^) is adopted here, in contrast to Ref. ^ In a proper analysis, Eqs. 
should be treated self-consistently, the problem which, to our knowledge, can be solved only 
numerically. Here, we are interested in a particular question of calculating the damping 
of QP hole states in an analytical way based on the well-established numerical results of 
SCBA. 

Actually, spectral characteristics of a hole propagating in an AFM background at low level 
of doping have been investigated in many workslll. Those results led to the consensus that 
the hole spectrum involves a narrow QP band of coherent states at low energies and a broad 
continuum of incoherent states. The corresponding spectral function is then represented as 
p^{u) = p^°\u) + p^^uj) with 

p'^''{uj) = Z^5{u; + p-E^). (6) 

The QP weight and the bandwidth W are estimated to be — J/t = 2" and W ~ 2 J. 
The QP dispersion in the vicinity of its minima kj = (±7r/2, ±7r/2) can be expanded as 
-f'ki+k' — Ei^i + ^|p/2^|| + A;^/2m_L.lii Here, fcy and k'j_ are the component of k — kj in the 
(1,-1) and (1,1) directions in the Brillouin zone (BZ) for k^ = (7r/2,7r/2). For instance, 
the anisotropy factor a = m\\/m±_ is calculated to be a = 6 for J = O.St.i This anisotropy 
can be absorbed by the following transformation k' {ak'p k'j_), which does not change our 
final results. Hence, we further consider the case my = m± = m(~ J~^).i^ 

From the above results one may expect that the filling of QP states leads to a four- 
pocket FS. Some arguments have also been given in Refs. ^ and |T2| that the fraction of 
BZ covered by these pockets at T = is equal to the hole concentration 6. This leads to 
the following estimations for the Fermi momentum kp = V^vr^ and the chemical potential 
/i = Ek, + k%/2m. 

The nearly structureless incoherent part pj^'^(u;) is distributed predominantly above the 
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QP band and can be approximated astj 

f^^\u') = {i/2V)e{\J\-j)e{2v-J), (7) 

where 2T ^ 2zt and uj' is measured from the middle of the QP band. In Eq.(|^ the neg- 
ative energy cutoff Uc — —J — 2r(l — Zk)5 is imphed. It is provided by the sum rule 
A^^^ / dujniuj) p]^{uj) = 6 taken at T = with the value of the chemical potential fi de- 
fined above. Numerical analysisiEl of the model allows us to conclude that Uc does not 
depend on T (for T <^ J). 

In the above formulated scheme (H)-® the QP damping is due to scattering by spin- 
waves. For QP states near FS renormalization of low-lying long-wavelength spin excitations 
[uj Ep = k'^p/2m, q ^ 2kp -C 1) is of crucial importance. This renormalization is due 
to the coupling of spin-waves to "particle-hole" pair excitations and is described by the 
polarization operator (|^) which contains three contributions. The first part n'^"'^(q, tu) is due 
to the transitions within the narrow QP band, when both pk(^^i) and Pk-q(^^2) in Eq-® 
are replaced by p™*^. The remaining two terms n'^"'(q, cu) and n^"'(q, u;) are provided by the 
coherent-incoherent and incoherent-incoherent transitions. 

First, considering I{^~^{(\^uj) we come to the following expression for small |q| <^ 1 

n:^(q, ^) = ^ E ^k,M„(k, + k', q)M^(k, + k', q) 

i,k' 

[^(gk'-q) - ^(gkO] .g^ 
UJ- Si,! + £k'-q + 

where Sk' = {k'"^ /2m ~ Sp) is the hole energy referred to the Fermi level Sp = kp/2'm and 
the summation over i is due to the presence of four equivalent minima. Since for small 
momentum |q| ^ 1 the vertex function Mi 2(k, q) is proportional to ^/g, we make an 
approximation Mi^2(ki + k', q) ^ ±M(kj, q), where 

M(k, q) = 2-5/4g-i/2 (g^ gin k, + qy sin ky), (9) 

to keep the leading contributions in Eq.(^. This leads to the following relations between the 
elements of the polarization operator: n5"]^(q, cu) = n22^(q, c<j) = —Yi'{2 oj) = —Il2i{q,u!) = 



nq(a;). We also note that the summation over i in Eq.(P) introduces an effective interaction 
|M(k„q)|2 = q/y/2. Then, for T = we obtain the following expressions for the real 
and imaginary parts of 

ReU,{io) = C{-q + sgn[r]{q,uj)][u{q,uj)]'/' 

+sgn[r/(g,-cu)][z/(g,-cu)]i/2|^ (^q^ 

Imn,(^) = C{[-z/(g,^)]i/2 - [-u{q, -u)]'/'}, (11) 

where rj{q,uj) = muj/q + q/2, v{q,uj) = 7f{q,uj) — {kpY and C = A\^mt'^ Z'^ / it ~ A\/2J/tt. 
The step-G-functions insuring the positivity of the arguments of the square roots are implied 
in Eqs.(|lOD and (0). 

For further purposes we fix also the asymptotic, u ^ 0, behavior of Ilq{uj) ioi q < 2kF'- 

Imn,(^) = 0, ReUqiu) = (12) 

for uj/q > kp/m, while for uj/q < kp/m one has 

—2Cmuj 

Imn„(cj) = , ReUqiiu) = -Cq. (13) 

\2kFy - g2 



The limiting case ([T^) is important in calculating of the renormalized spin-wave velocity, 
which will be shortly discussed below. The limit ( |T3|) corresponds to the region of the 
spin fluctuation spectrum, generated by "particle-hole" excitations, which produces finite 
damping of quasiholes near FS. In this respect, we note the conventional linear cj-dependence 
of lmllg{uj) in Eq. ([T3|) in contrast to a marginal (g- and cj-independent, at T = 0) form 
of Imng(ti;) in Ref. ^ (see Eq.(19) there). The reason of this difference can be explained as 
follows. Due to the incorrectly defined limits of integration in Eq.(14) of Ref. |^ the authors 
lost part of the polarization operator (which is presented by the second term in the curly 
brackets in Eq. ( [Tl| ) in this paper), which led to wrong subsequent approximations. 

Considering the remaining contributions n^"^(q, cu) and n^Jj(q, cu) we note that each of 
them is characterized by a threshold energy A for creating a "particle-hole" pair excitation. 
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Namely, A = ep ioi processes involved in n^"^(q, cu) and A = 2J for ]T^p{q,u!). Therefore, 
Imll™! = 0, where 11™^ = Il^'^ + nj^Jj, for frequencies u < ep are interested in. Evaluation 
of Rell™^ requires the summation over all virtual processes, which gives a finite estimate for 
these quantities even as u; — > 0. To the lowest order in q and 6 for Rell™! we obtain 
Renfi22i(l^^) — ^RsII'/f 2i(Q'? ^) — —Aq6 where A is positive and can be estimated as 
A ~ t/V2{\n{zt/J) + z^{l - Z k)[\n{2 J /zt6) + 1]}. 

A position of the pole in the spin-wave GF in the long-wavelength limit is now deter- 
mined as 



c^q = ujjl - 2[Aq6 - ReUg{u)]/ujg + 0{6'). (14) 



Since the unrenormalized spin-wave velocity u = \plJ is much larger than the Fermi velocity 
vp = \fn5J one has to take in Eq. ([T^ the limit (|^) which gives ReY\.q{ujq) ~ U0q5. Thus 



the renormalized spin- wave velocity now reads u{6) = uJl — 2{A/u — 1)5. For the actual 



values of 5 it holds A/u ^ 1 and hence, one obtains, in accordance with Refs. ^ and a 
spin-wave softening due to the presence of the incoherent part in the hole spectrum. 

The above estimation for u is valid up to the critical hole concentration 5c which is 
defined as u{5c) = vp{= ^/^l5^J). In particular, for J/t = 0.3 we estimated 5c — 0.04. 
For higher concentrations 5 > 5c, by taking the corresponding limit (^) one can see from 
Eq. ( |I4D that the pole uj^ becomes purely imaginary. So, the long-wavelength magnons, with 
q ^ 2kp lose their identity and can not be now detached from the incoherent part of the 
spectrum produced by pair excitations. Disappearance of the long-wavelength magnons 
due to dilution of the AFM state with holes was connected by several authors0 with the 
occurrence of a phase transition into a disordered magnetic phase. The applicability of the 
spin-polaron model in the disordered phase will be discussed bellow. 

Let us consider the effective spectral function (^. By using the approximated vertex 
function (j^) for the momentum near FS, k' kp (k' = k — kj), we obtain 

Xk.+k',qM - -l/7r|M(ki,q)|2lmD(q,a;), (15) 
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where M(ki,q) = 2~^/^g^/^qkj and D(q,uj) = Du — 2Di2 — D22- Taking into account the 
relation between the elements of the polarization operator n(q, cu), that has been obtained 
above, we write 

Im£)(q,a;) = AcullmU.ico) /\V^ico)\^ . (16) 

For the actual region of the uj- and g-variables, defined as u/q < vp, where the part 
of the spin fluctuation spectrum responsible for the quasihole damping is located, one has 
|r'q(u;)P = [cu^ + c^gY + [2uj^lmllg{uj)]'^ with c = u^2{A + C)/u- 1 > vp- This strong 
inequality allows us to take the static limit, u —>■ 0, for I^q(ci^) in Eq. (|l^) that results in 

Xk,+k',q(^) ^ -l/7rV2(M/c2)^Imn,(a;)(qk,)2g~i. (17) 

Inserting (|17D into @ one obtains for the imaginary part of the hole self-energy 

e 

lmS(k',£:) oc JJ dq cos^6d6jdujlmIlq{uj)S{e — to — ^k'-q), 



where cos 6' = qkj, k' = k — kj, £ is the hole energy referred to the Fermi level and Imnq(u;) 
is defined in Eq. (|l^). 

Like in the conventional considerations of the 2-dimensional (2D) FlJiI the major con- 
tribution to ImS is given by scattering processes with the momentum transfer q almost 
parallel to k'. These processes result in the following dependence for ImS familiar for the 
2D FL:0 

ImS(k',£) oc f{\k'ki\){eyeF) ln{e/eF). (18) 

Here, a k'-dependence of ImS is due to the anisotropy of the vertex (^ and is given by 
/(|k'kj|) which is a positively defined smooth function of its variable. 

Up to now we have considered the scattering processes retaining a quasihole in the vicinity 
of the same hole-pocket. There exist, however, processes in which the hole scatters from 
a given hole-pocket to the opposite or neighboring one, with momentum transfer q ~ Q 
(Q is the AFM wave vector) and q ~ Q' = (vr,0), respectively. Since the symmetry of the 
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problem provides the equivalence of g ^ 1 and g' = |Q — q| ^ 1, the first process gives just 
an extra factor 2 in ImS. Further, the vertex function M(kj, q) falls much faster at q ~ Q' 
than at q ~ (or Q) and the second kind of processes gives higher order corrections in uj. 
So, the conventional 2D FL behavior is expected for quasiholes at low doping. 

Being originally formulated for a state with an AFM ordered spin subsystem, the spin- 
polaron model requires some justification if one tries to extend it to a disordered phase, i.e. 
either to 5 > 5c at T = or T > 0. Actually,^ the hole spectrum is weakly affected by 
the absence of the long-range order, provided that the AFM correlations with the radius 
^ ^ Rp survive {Rp is the size of the spin-polaron associated with a hole). It means that 
hole propagation over the same sublattice dominates and the four-pocket FS survives as 
well. 

Connecting the magnetic phase transition aX 6 = 6c with the disappearance of long- 
wavelength magnons with q ^ 2kF, we did not find, however, any abrupt change in that 
low-lying part of spin fiuctuation spectrum which is responsible for the quasihole damping. 
Therefore, one may expect that not only the QP dispersion relation but also the character 



of the quasihole damping (TTBI) do not change for S slightly above 5c. 

This picture breaks down with further dilution of the magnetic subsystem, when the 
magnetic correlation length becomes comparable with the size of the spin-polaron. In this 
case, the nearest -neighbor hole hopping becomes dominant and a transition to a large FS 
takes place. However, this regime is beyond the scope of the present consideration. 

Considering a possible effect of finite T (low enough to provide ^ ^ -Rp ) we point out 
the existence of a characteristic temperature Td{S) above which one may expect different 
behavior of a quasihole subsystem as compared to the low temperature case, T -C Td{6). 
Actually, the Fermi-ensemble of quasiholes goes over into the strongly nondegenarate regime 
when the temperature Ta{S) ^ ep ^ 1.5J6 is reached (for instance, ~ lOOK at 5 = 0.05 
for J = 1500K). That is a result of the strong renormalization of the chemical potential n 
with T, which is naturally inherent in a fermion system of low density. Really, our analytical 
estimations, as well as numerical calculations,^ show that /i crosses the bottom of the QP 



band at T ~ and for T > Td lies in the low energy incoherent part of the hole spectrum. 
This results in a dramatic change in the momentum distribution function N{k): the four- 
pocket structure existing at T < is almost washed out at T > The onset of this 
strongly nondegenerate regime for hole carriers should manifest itself in a strong change 
of the thermodynamic and transport properties of the system, the problem which requires 
further theoretical and experimental studies. 

In summary, we have investigated the quasihole damping in the low doping regime, 5 <^ 1, 
of the 2D t — J model. The self-energy parts for the hole and magnon GF are derived within 
the self-consistent Born approximation. Based on the well established result Ji3 for the spec- 
tral density function of a hole moving in the AFM background, we first have calculated 
renormalization of spin-wave excitations due to the presence of holes. With increasing hole 
concentration 6, softening of the long- wavelength spin-waves followed by their overdamping 
aX 6 > 5c(~ 0.04) has been obtained. The renormalized spectrum of spin excitations was 
incorporated to calculate the imaginary part ImS(k, e) of the hole self-energy. It has been 
shown that ImS, as £ 0, possesses the form ( ]TB| ) characteristic of the conventional 2D 
PL. 
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